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Abstract. This article presents the study of the light-curves extracted from the MACHO database of a sample 
of stars observed by the Infrared Space Observatory in the Small Magellanic Cloud. These stars belong to the 
ISO-Mini-Survey catalogue of the Magellanic Clouds (ISO-MCMS, Loup et al. in preparation). Most of them are 
in the asymptotic giant branch (AGB) and supergiant phases. The dominant period and amplitude of pulsation 
have been derived and the stars have been classified as Mira or Semi-Regular pulsators. Furthermore, the cross- 
identification with near-infrared DENIS and 2MASS magnitudes available within the ISO-MCMS allowed us: (i) 
to investigate the properties of these stars in the combined near- and mid-infrared colour-magnitude diagrams, 
(ii) to derive the bolometric magnitude by integrating the spectral energy distribution and (iii) to estimate the 
mass-loss rate. The stars have been divided into carbon- (C-) and oxygen-rich (O-rich) using the (J — Ks, Ks) 
colour-magnitude diagram and their period and amplitude distributions have been compared. C-rich AGB stars 
have a sharp peak in their period distribution at about 250 days and have on average a larger amplitude than 
O-rich AGB stars. This effect, not previously detected from the study of similar stars in the Large Magellanic 
Cloud and in the Baade's window, might be closely related to the metallicity of the environment in which the 
stars have formed. 
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1. Introduction 

The stellar content of the Small Magellanic Cloud 
(SMC) can be studied in great detail because of the 
low galactic extinction and a known distance m— M= 
18.99±0.03(formal)±0.08(systematic), Cioni et al. il2UUU|) . 
However, it is well known, but not yet quantified, that the 
SMC is extended along the line of sight much more than 
its companion galaxy, the Large Magellanic Cloud (LMC). 
Together with our own Galaxy, they constitute an inter- 
acting system (e.g. Gardiner & Noguchi 119961 Moore & 
Davis 119931 Heller & Rohlfs IT^M|l and a metallicity scale 
environment in which the Galaxy is the most metal rich 
and the SMC the most metal poor of the three. Fields cov- 
ering almost the whole of the two Magellanic Clouds and a 
considerable number of fields in the Galaxy have been re- 
peatedly observed by the micro-lensing projects: OGLE, 
MACHO and EROS. This is providing a wealth of data 
to identify and study variable stars in every evolutionary 
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phase. Partial or all sky survey data in the near-infrared 
(IR) provide essential counterparts for the specific ana- 
lysis of the reddest objects: red (RGB) and asymptotic 
giant branch (AGB) stars. Most AGB stars are pulsating 
with an amplitude of up to a few magnitudes in the optical 
wave bands and somewhat less in the near-infrared bands. 
Long Period Variables (LPVs) have a period of at least 
100 or 150 days and longer. Many are also multi-periodic 
objects. More than one dependence between period and 
luminosity has been detected in the LMC (Wood et al. 

Cioni et al. l2TiuTl Noda et al. l^uTH Lebzelter et al. 
2002|. Each period-luminosity relation (PLR) is proba- 
bly associated with a different mode of stellar pulsation. 
Because most of these stars are probably experiencing the 
thermal pulsing AGB phase (TP- AGB) they are losing 
mass at an average rate of 10 _6 MQyr _1 and they have 
highly extended atmospheres. The pulsation as the driving 
mechanism of the stellar wind is a fundamental element 
in the study of the mass-loss rate of AGB stars. They 
are surrounded by an initially optically thin, and later op- 
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tically thick, circumstellar envelope and thus sometimes 
become detectable only in the infrared wave bands. Many 
obscured AGB stars were discovered by the IRAS satel- 
lite. More recently the Midcourse Space Experiment (MSX 
- Price et al. I2001[) observed selected high density regions 
in the Magellanic Clouds. Improved source position were 
derived. Egan et al. I|2001fl present the cross-identification 
between MSX and 2MASS data in the LMC. However, 
the Infrared Space Observatory (ISO) obtained more sen- 
sitive data in some fields in the Magellanic Clouds and the 
Galaxy that increased the number of obscured AGB stars 
and allowed us to learn in more detail about their dust 
properties. 

In order to study the evolved AGB population it is 
thus important to utilize observations that cover the spec- 
tral range from the optical (/ band) to the mid-IR. It 
is then possible to derive fundamental quantities such as 
the bolometric magnitude and the mass-loss rate of the 
stars. In addition, from the period, the amplitude and 
the regularity of the variation, it is possible to charac- 
terise the stars as Mira or Semi-Regular (SR) pulsators 
and to indicate, from the PLR, the mass of the stellar 
progenitor. In this article the properties of a sample of 
AGB stars in the SMC extracted from the ISO-Mini- 
Survey catalogue of the Magellanic Clouds (ISO-MCMS, 
Loup et al. in preparation), which already contains the 
cross-identification with the near-IR catalogues DENIS 
and 2MASS, are discussed. This is done by additionally 
analysing the light-curves extracted from the MACHO 
database. A similar work was performed in the LMC bar- 
west field (Cioni et al. 12001(1 and will be further extended 
to other fields were ISO measurements are now available. 
A detailed comparison of more or less obscured AGB va- 
riables between the two Clouds and the Baade's window 
NGC6522 will be given in another paper (Glass, Schultheis 
& Cioni, in preparation). 

The data and the cross-identification with the light- 
curves are described in Sect. 2. In Sect. 3 the analysis of 
the period and amplitude of variability and the determi- 
nation of the bolometric magnitude is presented. A dis- 
cussion of the different type and variability class of the 
sources in the PLRs and in the near- and mid-IR colour- 
magnitude diagrams is given in Sect. 4. Sect. 5 concludes 
the article. Appendix A discusses the re-analysis of the 
LPVs studied by Feast et al. (|1989[l and B compares the 
bolometric correction obtained in this work with the rela- 
tions by Alvarez et al. |20Q01l and Montegriffo et al. lfT33%|) . 

2. Data 

The sample of data in the SMC, analysed here, is based on 
the cross-identification between the ISO-MCMS (Loup et 
al. in preparation), and the MACHO light-curves publicly 
available 1 . 

The ISO mini-survey in the direction of the SMC co- 
vers an area of 0.28 square degrees. Observations were per- 



formed with the infrared camera ISOCAM on board the 
ISO satellite in the LW2 (6.61/mi) and LW10 (11.18/xm) 
wavebands. The catalogue lists 1333 sources brighter than 
about 12.4 and 11.3 mag in LW2 and LW10, respec- 
tively. The coordinates of the ISO sources are better than 
4" with a faint tail up to 6". The catalogue also lists 
the cross-identification with the near-infrared catalogues 
2MASS (2nd Incremental Release PSC 2 ) and the DENIS 
(Epchtein et al. [1997) catalogue towards the Magellanic 
Clouds (DCMC - Cioni et al. I2000a|l . and a few smaller 
optical catalogues of carbon stars, HII regions, PNe and 
emission-line stars. According to the association likeli- 
hood criteria defined by Loup et al. (in preparation) 
48% of the sources have no, or an unlikely, association 
with DENIS/2MASS; 13% have a questionable associa- 
tion; 25% have a confident association and the remaining 
14% are "known" sources with a DENIS/2MASS and op- 
tical counterpart. It is perhaps surprising that there are 
many ISO sources without a near-infrared counterpart. 
However, most of these sources are among the faintest ISO 
detections in both ISO wavebands. The histogram of their 
magnitudes and colour peak at: LW2 = 12, LW10 =11 
and (LW2 - LW 10) = 1.5. Moreover, the ISO-MCMS 
includes sources well detected in one band only. Among 
sources detected only in two bands about 86% have a 
near-IR counterpart. The ISO photometric magnitudes 
are part of the final version of the catalogue released to 
the Co-Is and associates on March 2002. 

The MACHO (Alcock et al., 2000) light -curves were 
obtained from repeated measurements in two, blue and 
red, broad pass-bands from the 50 inch telescope at the 
Mount Stromlo Observatory in Australia. Observations of 
fields covering most of the Magellanic Clouds were per- 
formed from 1992 to 2000. The light-curves are easily ac- 
cessible via the web interface giving in input the right 
ascension (RA) and the declination (DEC) of each source. 

Fig. n shows the regions observed by ISO superim- 
posed on the MACHO fields. The area covered by ISO is 
entirely contained in MACHO fields number 211 and 212; 
a few sources are also present in field number 207. These 
MACHO fields were observed for the whole period of eight 
years while the observations of fields in the external parts 
of the SMC cover just five years. The OGLE light curves 
cover a period of about five years (observations started in 
1997). Because of the nature of LPVs, sources with periods 
from 100 to 1000 days, the MACHO light-curve database 
was preferred. 

2.1. Match between ISO-MCMS and MACHO 

Each ISO source with a confident DENIS/2MASS as- 
sociation has been searched for a light-curve in the 
MACHO database using, as input, the 2MASS coordi- 
nates. Variables were confidently identified within a radius 
of 1'. Fig. shows that most of the sources were actually 
associated within 1". In a very few cases the MACHO 
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0.55 
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Fig. 1. Location of the fields observed by ISO (continuous 
line) and MACHO (dashed line) in the SMC. The ISO 
fields lie in the bar of the galaxy and are each 16' x 16' 
in size. Both right ascension (RA) and declination (DEC) 
are in degrees. 



Fig. 2. Histogram of the absolute distance between the 
2MASS counterpart of an ISO-MCMS source and the cor- 
responding MACHO source. 



counterpart, a source with a typical LPV light-curve, was 
not the closest to the 2MASS source - but still within no 
more than 4". We consider a light-curve to be that of a 
LPV when a more or less periodic behaviour is detected 
by eye. It is expected that most of the sources detected in 
the near- and mid-IR are associated to AGB stars. This is 
a consequence of the mass-loss process occurring in AGB 
stars that builds up a circumstellar envelope observable in 
the mid-infrared. 

Among the 502 sources with a confident 
DENIS/2MASS counterpart, 458 have been found 
to have a MACHO light-curve. Table electronically 
available, lists the MACHO coordinates at the epoch 
J2000 [Column 1 & 2), the MACHO identifier field, tile 



and number [Column 3, 4 & 5), the distance in arcsec 
between the 2MASS and the MACHO source (Column 6) 
and the ISO-MCMS identifier [Column 7). The first 10 
lines are given as an example (Table 

The sources listed in Table |21 are those for which we 
could not find a convincing MACHO counterpart. The 
ISO-MCMS name is given in Column 1, the distance in 
arcsec to the first possible MACHO counterpart is given 
in Column 2 and in some cases the reason of the miss- 
ing cross-identification is written in Column 3. Too few 
data points means that there are not enough points to 
characterise the shape of the light-curve. This would be 
useful to recognize a MACHO counterpart within the ISO 
astrometric precision. 
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Table 2. Table of ISO-MCMS sources with a confident 
DENIS/2MASS counterpart but without an associated 
MACHO light-curve. 



ISO - MCMS MACHO counterpart ? 



J004850.9 




TO -l A AA 

7314U2 


a > 17 


TAA A A A A r* 

J004902.7 




TO A 1 A 1 

732141 


J ^ C" 1 O J J- "J- 

d > 5 , only 3 data points 


taa a a a n a 

J004907.2 




TO 1 O O f 

731335 


J ^ 1 a// J- X* l j_ j_ 

a > 12 , too lew data points 


J004907.8 




730920 


d > 14" , too few data points 


J004908.2 




731415 


a > 14 


J004909.5 




TO 1 1 O O 

731133 


d > 7 


TAA ,1 A 1 A f 

J004910.5 




TO "1 Of A 

731859 


d > 3 , too lew data points 


TAA A no A A 

J004934.4 




to i a nn 

731409 


d > 9 


Taa AC\ AC\ O 

J004940.8 




TO 1 oco 

731352 


d > 15 , too lew data points 


JC\C\ AC\ AC\ f\ 

J004940.9 




TO 1 A 1 O 


d > 8 


J004959.0 




731058 


d > 7" 


J005001.7 




TO 1 o a n 

731349 


d > 9 


J005034.2 




TO 1QC/1 

731354 


7 ^ i r// 

d > 15 


taa i — i r\ o i 

J005108.1 




TO 1 O A O 

731343 


d > 5 


taa r -i -i t r 

J005117.5 




to i^cn 

731659 


d > 4 , too few data points 


TAA 1 1 A 1 1 

J005121.1 




731416 


d > 4" 


TAA 1 1 A A T 

J005129.7 




TO 1 A A A 

731044 


d > 8" 


TAA 1 1 O 1 A 

J005131.4 




TO AAAO 

732008 


d > 5", too few data points 


TAAI — 1 1 a A 

J005140.2 




TO 1 O O O 

731333 




J005202.5 




TO 1 O O A 

731339 


J ^ all 

a > 6 


J005207.0 




TO 1 /in7 

731407 


d > 7 


TAA f A 1 1 f 

J005211.5 




731351 


d > 16 


TAAfAl T O 

J005217.8 




TO 1 OAT 

731327 


d > 4" 


J005218.9 




/ 30921 


d > 4" , too few data points 


TAACAOA O 

J005230.8 




TO AA f A 

730259 


d > 8 , too lew data points 


J005239.9 




Tonn i o 

730912 




taa ro^C T 

J005245.7 




730425 


7 . r /l 

a > 6 


TAAfOAA A 

J005302.9 




TO AT 1 /"' 

730746 


d > 6 


J005306.4 




7onc o o 

730633 


d > 10 


TAACOAA A 

J005309.0 




TO A A f\ A 

730404 


d > 9" 


TAACOAA A 

J005309.9 




TO 1 O A A 

731324 


J ^ o" j_ X* 1 j_ • j 

d > 3 , too lew data points 


7r\r\ coin o 

J005310.3 




TOriTOO 

730722 


a > o , too tew data points 


TAA I"" ,1 A |— If 

J005405.7 




TO A 1 f A 

730459 


d > 3 , too lew data points 


J005417.7 




730504 


d > 3", poor sampling 


J005418.1 




730523 


d > 7" 


J005518.6 




731420 


d > 2", too few data points 


J005549.0 




730257 


d > 4", too few data points 


J005646.4 




730519 


d > 12" 


J005648.0 




730404 


d > 5" 


J005650.1 




731040 


d > 2", too few data points 


J005650.9 




730420 


d > 7" 


J005749.0 




730522 


d > 11", too few data points 


J005822.0 




730734 


d > 4" 


J005833.1 




730528 


d > 16" 



3. Period analysis 

The Schwarzenberg-Czerny method (1996) has been used 
to detect periodicities in the light-curves of the selected 
sources, independently for the measurements in the blue 
and red bandwidth of MACHO. The algorithm was exe- 
cuted thirty times searching for periods in a range from 10 
to 100, 200, 3000 days with a step of 0.005 days. The re- 
sult of the program is the period of the variation and a few 



other parameters among which the quality of the variabil- 
ity "q" . We identified as the main period the value corre- 
sponding to the highest value of "q" among the thirty dif- 
ferent determinations. We noticed that at the first passage 
through the periodic detection algorithm the dominant 
period did not give, in some cases, a sufficiently defined 
phase-curve (distribution of the measurements folded on 
a given period). By restricting or enlarging the range of 
searching periods it was always possible to optimize the 
phase curve to a certain period. Therefore we decided to 
execute the algorithm more times by changing the pe- 
riod searching window. The reliability of the assigned pe- 
riod was then inspected by looking at the phase curve 
of each source. Sources of regular variability and with a 
large amplitude of variability have a clean phase plot with 
a large S/N ratio that allows the variation of the posi- 
tion of the maximum or minimum brightness and/or the 
presence of bumps to be distinguished (sometimes identi- 
fied as the secondary period). In other cases the S/N was 
not good enough to detect these secondary variations but 
the phase-curve was clear enough to confirm the main pe- 
riod. Then we introduced the following criteria to flag the 
sources in each band. Flag= for sources that do not have 
data points in a given MACHO band. Flag= 5 for sources 
with q < 15 or a period that coincides with the extremes 
of the searching window in all 30 cases (i.e. sources with 
a period of 10 or 100, 200,... 2900, 3000 days). Their light 
curve is judged to be of poor quality. Flag= 1 and Flag= 2 
for sources with q > 15 and for which the assigned period 
(P) differs not more than 10 days from the mean of all the 
measurements (among the thirty possible) with q > 15 
in the blue and red band width, respectively. With these 
criteria we are confident that a single periodicity is the do- 
minant source of variation in the light-curve. Flag= 6 for 
sources with q > 15 and a difference from the same mean 
larger than 10 days. These sources are likely to have more 
than one significant periodicity. In this work we are not 
aiming at a detail analysis of the secondary periods which 
will be discussed, together with the comparison of LPVs 
with an ISO counterpart between the Magellanic Clouds 
and the Galaxy , in Glass, Schultheis & Cioni (in prepa- 
ration). Aliases were identified from the diagram Log(P) 
versus K$ as those periods creating clear vertical patterns. 
These correspond to periods equal to 29 days and in the 
range from 340 to 390 days, included. These sources were 
assigned Flag= 9. Examples of MACHO light-curves for 
LPVs can be found in Wood et al. <{l"39"9"|) . Table El sum- 
marises the detected periodicities. It lists the description 
of the Flag ( Column 1) , the value of the flag and the num- 
ber of sources corresponding to a given flag in the red 
(Column 2 & 3) and blue bandwidths (Column 4 & 5). 

Table 01 provides the magnitudes from the different 
databases, the bolometric magnitude (see next Section) 
and the parameters that describe the dominant pulsation 
of all the sources in the sample with an extracted light- 
curve (458). This table is also made electronically avail- 
able and the first 10 lines are given as an example (Table 
EI). The content of the table is the following: ISO-MCMS 
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Table 3. Summary of detected periodicities 



Description 

No curve 

Good 

Poor 

Multi — Periodic 
Alias 



Flag N Flag N 



RED 



BLUE 



8 



9 



96 2 91 

103 5 208 

91 6 91 

160 9 59 



name ( Column 1) , IJKs DCMC magnitudes ( Columns 2, 
3 & 4), JHK S 2MASS magnitudes {Columns 5, 6, & 7), 
LW2 and LW10 ISO magnitudes {Column 8 & 9), the ap- 
parent de-reddened bolometric magnitude {Column 10), 
period, amplitude, flag and quality parameter derived for 
the light-curve in the MACHO-red bandwidth {Column 
11, 12, 13 & 14) and in the MACHO-blue bandwidth 
{Column 15, 16, 17 & 18). 

Values of 99.99 for the magnitudes and of —99 for pe- 
riod and amplitude indicate missing data. Dclmotte ct al. 
12002 have shown that Jdcmc = J-ZMASS ~ (0.11 ± 0.06) 
and K sDC MC = K s2M ass - (0.14 ± 0.05) but we did 
not apply this correction to Table 4. Amplitudes (A) 
are the difference between the minimum and the maxi- 
mum value of the MACHO intensity and are expressed in 
mag x 100. In a handful of sources this calculation might 
be biased by spurious data points with a magnitude ex- 
cursion compared to the sinusoidal behaviour of the light- 
curve. Perhaps these points are real and represent a flare. 
The period is in days. 

3.1. Bolometric Magnitude 

With the availability of magnitudes from the / band to 
the mid-IR band LW10, it is possible to determine quite 
precisely the bolometric magnitude of the sources in the 
sample by integrating under the spectral energy distri- 
bution (SED). This integration has been done by fitting 
a spline to the SED profile in the plane {v, F v ), extra- 
polating linearly to zero flux at zero frequency, and also 
extrapolating linearly at the other extreme of the SED 
through the two near most points. The SED of all the 
sources detected by DENIS and ISO were inspected one- 
by-one to exclude those which are either too red or too 
blue to describe satisfactorily their SED with the available 
measurements. 

The adopted extrapolations represent the SED very 
well, almost independently from the circumstellar opacity 
for all M spectral subtypes (Loup et al. in preparation). 
That article gives a detailed determination of the multi- 
band bolometric correction for AGB stars by combining 
data from the major survey instruments: DENIS, 2MASS, 
ESO, IRAS, ISOCAM and MSX. The bolometric correc- 
tion is obtained from the different combination of pass- 
bands and a theoretical model for the SED. This is com- 



pared with the correction obtained from the integration 
under the SED using a spline or a linear fit through the 
flux points. A fairly well sampled SED, that detects most 
of the energy emitted by the objects in a set of characte- 
ristic bands, is a necessary condition to derive good bolo- 
metric luminosity {Mbol) for stars of different M spectral 
subclass and C stars. 

The best wavelength coverage for the data in the sam- 
ple studied here is obtained by combining DENIS and ISO 
magnitudes. Though 2MASS K$ magnitudes exist for 455 
sources, DCMC-JiCs magnitudes exist for 366 sources 
and DCMC-7J mag exist for about 450 sources. On one 
hand it would seem natural to combine I J DENIS mag- 
nitudes with 2MASS-ifs magnitudes and the ISO magni- 
tudes, for the same sources, to derive the bolometric mag- 
nitude. On the other hand the DCMC~IJKs magnitudes 
were obtained by simultaneous observations and this is a 
unique constraint to the SED profile. 

Note that the amplitude of variation in the mid-IR 
bands can be as high as 1 mag (e.g. van Loon et al. 1998 
van Loon et al-EDOl), this effect may influence the shape 
of the SED during the variability cycle. The constraint of 
the I band to the shape of the SED at higher frequencies 
is very important for optically thin AGB stars because the 
peak of the SED is at about K$- In Loup et al. {in prepa- 
ration) we show that the linear extrapolation to zero flux 
at the longer frequencies starting from the flux in the J 
band subtends a larger area than the extrapolation start- 
ing from the flux in the / band. This effect overestimates 
the bolometric magnitude. 

Prior to the calculation of the flux emission in the dif- 
ferent bands and the integration under the SED we de- 
reddened the data adopting the extiction law by Glass 
For the DENIS pass-bands [A v : Ai : Aj : A Ks = 
1 : 0.592 : 0.256 : 0.089], and we assumed zero absorp- 
tion in the ISO wave bands. We used Ry = 3.1 and 
E{B — V) = 0.065, the latter is the average of the mea- 
surements, in the SMC, discussed in Westerlund ((T997 ) . 
The zero-point of the bolometric scale is computed with: 
M BOL ,Q = 4-74 and / iot , = 1.371 x lO" 8 erg cm' 2 s" 1 . 
A comparison of our bolometric magnitudes with the re- 
lations obtained by Alvarez et al. 12110.(3) and Montegriffo 
et al. (1998) is given in Appendix B. 

4. Discussion 

4.1. Period-luminosity relation (s) 

Fig-Et shows the relation between the logarithm of the pe- 
riod (days) and Kg. Variable stars are distributed on four 
different parallel relations. These relations were first dis- 
covered by Wood et al. ( 1999) from the analysis of a sam- 
ple of MACHO stars in the LMC. They were subsequently 
found by Cioni et al. I|2001|l using EROS2 (Lasserre et 
al. EOOO) light-curves, by Noda et al. lj2"o"o"2")> using the 
MOA database and by Lebzelter et al. 1|2002[1 using the 
AGAPEROS data. These latter measurements refer also 
to samples of stars in the LMC. A preliminary compar- 
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Table 4. Extract of the table of the properties of LPVs with ISO counterpart in the SMC 



ISO - MCMS 


I 


J 


Ks 


J 


H 


Ks 


LW2 


LW10 


niboio 


P 


A 


f 


q 


P 


A 


f 


q 








DCMC 




2MASS 


ISO 




M. 


- RED 


M. 


-BLUE 


/004844 4 


— 732119 


14.69 


12.91 


99.99 


13.09 


12.11 


11.80 


10.96 


9.84 


99.99 


203 


428 


1 


70 


374 


74 


9 


67 


J004844.5 


- 731754 


14.65 


12.82 


99.99 


13.18 


11.92 


11.14 


10.18 


9.84 


99.99 


343 


91 


9 


253 


342 


180 


9 


277 


J004844.7 


- 732018 


9.65 


8.57 


99.99 


8.74 


8.07 


7.99 


7.76 


8.06 


99.99 


2900 


401 


5 


20 


2900 


334 


5 


15 


J004849.4 


- 731627 


14.29 


13.99 


99.99 


14.19 


14.09 


13.77 


11.64 


10.91 


99.99 


1900 


95 


5 


46 


2063 


82 


6 


116 


J004849.8 


- 732002 


11.80 


10.31 


99.99 


10.48 


9.65 


9.42 


9.26 


9.24 


99.99 


29 


220 


9 


10 


821 


226 


5 


10 


J004850.7 


- 732138 


14.68 


13.67 


12.68 


13.80 


13.15 


13.05 


11.77 


99.99 


99.99 


368 


24 


9 


31 


1667 


42 


5 


11 


J004851.8 


- 732239 


10.90 


9.45 


8.41 


9.63 


8.83 


8.57 


8.04 


7.52 


11.21 


352 


456 


9 


46 


350 


301 


9 


7 


J004854.0 


- 731557 


13.62 


12.37 


99.99 


12.59 


11.85 


11.69 


11.58 


99.99 


99.99 


10 


19 


5 


68 


10 


21 


5 


19 


J004859.7 


- 731155 


13.04 


11.80 


99.99 


11.96 


11.17 


11.00 


10.74 


99.99 


99.99 


368 


184 


9 


21 


21 


29 


5 


4 


J004900.4 


- 732224 


14.16 


12.47 


99.99 


12.65 


11.70 


11.35 


10.74 


99.99 


99.99 


239 


67 


1 


137 


237 


99 


2 


207 



ison between the relations in the LMC and in the SMC 
has been shown by Cioni et al. (j2002a), concluding that 
sources in both galaxies occupy the same relations in the 
diagram (Ks, Log(P)). 

There are a few differences that depend on the way 
the two samples were selected and the light-curves ex- 
tracted: in the SMC there are sources with longer periods 
because the MACHO data cover a longer time range than 
the EROS2 data. Thus it is possible to identify longer pe- 
riodicities. Also there is a relative lack of SMC sources 
around the tip of the red giant branch (TRGB) because 
the LMC sources were selected from the DCMC and not 
among those sources detected at 7 and 12 /xm. It may be 
that AGB stars that are not yet in the TP-AGB phase 
or are in a very early phase of their TP-AGB evolution 
have a mid-IR emission below the ISO detection limit. 
The ISO survey was indeed designed to find mass-loss, 
not to detect sources without circumstellar extiction. 

The analysis of the light-curves of LMC sources de- 
tected by ISO has just started. The Kg magnitudes are 
always a one epoch measurement, and coupled with the 
amplitude of variation, that in the Ks band can be as 
high as 1.5 mag for OH/IR stars, prevent us from detec- 
ting a difference (if any) in the width of the PLRs due to 
metallicity. The same situation holds for stars in the Bulge 
(Alard et al. 2001), which have the additional disadvan- 
tage that the distance to the AGB stars is too uncertain 
to find all four PLRs. Furthermore, AGB stars in stellar 
clusters are too few to outline the presence of one and 
possibly two PLRs. There is evidence that in the elliptical 
galaxy NGC5128 one and maybe two PLRs exits for LPVs 
fReikuba 12002)1 . 

The atmosphere of AGB stars is dominated by 
molecules of TiO and VO if the star is O-rich (C/0< 1) or 
by molecules of CN, C2 and other carbonaceous molecules 
if the star is C-rich (C/0> 1). Stars of both types are eas- 
ily statistically distinguished in the (J — Ks, Ks) colour- 
magnitude diagram (Fig. because of the different 
strength of the molecular absorption bands at these par- 
ticular wavelengths (Cioni et al. 120011 Loup et al. 2002). 
The dividing line at (J — K s ) — 1.33 takes into account 



the shifts calculated by Delmotte et al. I|2002|l between 
2MASS and DCMC magnitudes and has been defined in 
Cioni & Habing (2002) for the whole AGB population of 
the SMC. It was also formerly checked by Cioni et al. 
(2000J via the cross-identification with the large sample 
of spectroscopically identified C-rich stars of Rebeirot et 
al. ( 1983). Filled symbols in Fig.|3]indicate C-rich objects. 
In Fig. [21d they define relation B and C, and dominate re- 
lation D (possibly offset from the oxygen-rich stars). 

4.2. Mass and Age 

According to the Vassiliadis and Wood (1993.) prescription 
for the evolution of low- and intermediate-mass stars car- 
bon stars are less massive than about 3Mq . They form af- 
ter a certain number of shell flashes when the atoms of car- 
bon overnumber the oxygen atoms. Hot bottom burning 
may prevent massive stars from becoming C-rich. Each 
track drawn in Fig. Oh starts at the first shell-flash cy- 
cle on the TP-AGB phase. Most of the C-rich stars in 
the SMC sample have a mass between 1 and 3Mq. The 
brightest stars (Mbol < —7) are probably supergiants 
(Wood et al. 1992). The most massive AGB stars with a 
mass of about 5Mq are about 0.1 Gyr old. The oldest 
detected AGB stars are about 9 Gyr old, but most of the 
sample is from 0.6 to 2 Gyr old. 

In the LMC (Cioni et al. I2001[) C-rich stars occupy 
only the brightest part of the PLRs, thus in a lower metal- 
licity environment, such as that of the SMC, C-rich stars 
form at lower masses, unless these stars are in the lumi- 
nosity dip that follows a thermal pulse or at the minimum 
of their light-curve. Menzies et al. (2002) recently showed 
that the top magnitude of the AGB stars in Ks in the 
dwarf galaxy Leo I ([Fe/H]=-2.0) are mostly carbon stars. 
They have J - K s > 1.1. The value of the J - K s that 
discriminates between O-rich and C-rich stars is a func- 
tion of metallicity. The lower the metallicity, the bluer the 
colour. Obscured AGB star candidates (J—Ks > 2) might 
also be of C type (50% level) and the observation at two 
different epochs indicates that at least one of these object 
might be a Mira variable (AA's = 0.87). 



M.-R. Cioni et al.: ISO LPVs in the SMC 



7 



A B C D . 

— o ° ° ° 

TRGB o 

I .... I .... I .... I ... . 



1.5 2 2.5 3 3.5 
Log(P) 




(c) 



-7 7 

-6 r 



-4 TRGB 




-3 -4 -5 -6 -7 



1.5 2 2.5 3 3.5 
Log(P) 



Fig. 3. Diagrams for the sources in the sample with a 
MACHO-RED light-curve and / = 1 or / = 6. (a) K s 
versus Log(P) relation. Letters ABCD identify the paral- 
lel sequences as defined in Wood ct al. (1999). Filled sym- 
bols indicate probable C-rich objects, (b) PLRs and over- 
plotted theoretical tracks by Vassiliadis & Wood (1993). 
The numbers indicate the initial mass of the stars in each 
track, (c) Amplitude distribution as a function of bolomet- 
ric magnitude, (d) PLRs for sources with an amplitude 
above the dashed line indicated in (c), q > 100 and / = 1. 
The continuous line is the Mira PLR fitted through the 
crosses corresponding to P < 420 days obtained by Feast 
et al. I|1989fl ; a distance modulus of m — M= 18.99 (Cioni 
et al. I2000fl has been applied. Symbols identify O-rich 
stars (open circles), C-rich stars (filled circles) and Feast 
et al. 119891 variables (crosses). 



4.3. Large Amplitude Variables 

Fig. shows the amplitude of the light-curves as a func- 
tion of Mbol for C-rich and O-rich objects. The sources 
with A * 100 > 110 and / = 1 have q > 100 and are dis- 
tributed in the PLR(s) as shown in Fig. 01. They occupy 
relation C which is the location of Mira variables in the 
LMC and the Galaxy. Over plotted is the sequence de- 
rived by Feast et al. (1989) and the single LPVs (crosses) 
as re-analysed in this work (Appendix A). This sets the 
minimum level of the quality parameter for Mira variables. 
This sequence was obtained including both O-rich and C— 
rich AGB stars with a period P < 420 days. A distance 
modulus of 18.99 has been used (Cioni et al. I2000|l . The 
single carbon at Log(P) = 2.1 and Mbol = —6 is of large 
amplitude because of flare-like events super imposed on 
a regular variation of smaller amplitude, thus this source 
should not be regarded as a Mira star. 




Fig. 4. Variable O- and C-rich AGB stars (a). The his- 
togram of the amplitude of variation (b) and that of the 
period distribution (c) are shown for sources of both types. 
The amplitude as a function of period is shown in (d). 



The scatter in the PLR, of the sources in this study, re- 
sults from by the unavailability of bolometric light-curves. 
However, it is consistent with the sigma associated with 
the relation derived by Feast et al. 119891 Estimating the 
magnitude of this effect would be valuable for studies that 
plan to use AGB stars as tracers of the star formation his- 
tory of systems as far as 250 kpc (Kucinskas et al. 2002). 
This preliminary investigation on the use of AGB stars 
concerns mainly early-type AGB stars which vary with 
amplitudes much below that of Mira-type stars. On the 
other hand Kucinskas et al. combined broad-band obser- 
vations from the literature. It should not be neglected that 
the amplitude of variability of AGB stars is different in 
each photometric band. Only simultaneous observations 
in different bands, as those that will be available with 
GAIA, provide a coherent SED from which a single-epoch 
bolometric magnitude can be derived. 

Finally we conclude that a large amplitude and a large 
value of q indicate a clear and regular light curve; thus a 
Mira star. There are some sources of large amplitude and 
with a large value of the quality parameter suspected of 
multi-periodicity Their main period does not correspond 
to relation C. A few Mira variables in sequence D have 
been found by Alard et al. H2001|l . On the contrary, this 
PLR could be populated by binary stars with an AGB 
companion (Wood et al. 1999). In what follows we mainly 
discuss the properties of the stars with a light-curve in 
the red band of MACHO and / = 1 or / = 6. 
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4.4. Mira & Semi-Regular Variables 

The distinction between Mira and Semi-Regular (SR) 
stars has become more and more subtle with the avai- 
lability of more sensitive long term observations and large 
samples of stars. The stellar light-curve can be as regular 
in SR as in Mira stars. Lebzelter et al. l|2002ll define a 
classification criterium for LPVs based only on the regu- 
larity and multi-periodicity of the light-curves. Most of 
the Mira stars fall into the class of regular variables but 
not all regular variables have large amplitudes of varia- 
tion. SR belong to a second class and have on average a 
smaller amplitude than the stars in the first class. Into 
the third class fall irregular variables, that do not show a 
clear periodicity. The Mira phase is believed to be just a 
moment during the AGB evolution and the star can thus 
become or have been SR. This scenario is supported by 
observations of very nearby variable stars over a period 
of about 90 years (i.e. Percy & Bagbv IT9991 Kiss et al. 
2000). However the evidence, e.g. from globular clusters, 
suggests that most SRs are in the pre-Mira stage. 

Mira variables lie on one PLR (Feast et al. I1989[) and 
SR occupy this and other PLRs according to their mode 
of pulsation (Wood et al. Cioni et al. IMTT1 Noda 

et al. l2"UU2l Lebzelter et al. |2"UU2"|) . If Mira variables are 
selected as discussed in Sect. 4.3 then they are clearly 
concentrated on the C PLR. This selection criterium is 
based on the amplitude of variation. However, the sources 
are the most regular as well (q > 100). 

Whitclock et al. I|1997|) . by analysing the light-curves 
of a dozen C-rich AGB stars in the Milky Way Galaxy, 
showed that the distinction between Mira and SR varia- 
bles is not necessarily very clear among carbon stars, and 
may not be useful as a distinction for C-rich stars as it 
is for O-rich stars. Fig. 0Ja shows the histogram of the 
distribution of the amplitude of variation separately for 
O-rich and C-rich stars. Perhaps there is a distinction 
between Mira, the pulsators at the tail of the amplitude 
distribution(s), and SR. Both O-rich and C-rich sources 
have a histogram with a dominant peak at about Ax 100 = 
40 and A x 100 = 80, respectively, and a tail to larger 
amplitudes. The fact that C-rich stars have on average 
larger amplitudes than O-rich stars is different from what 
was found in the LMC. There, both histograms had the 
same distribution (Fig. 4 in Cioni et al. (2001)). The same 
histogram but for the amplitudes derived from the BLUE 
light-curve shows that O-rich stars peak at about A * 
100 = 50 and C-rich stars at about A * 100 = 100. This 
confirms that the average amplitude is larger in the blue 
than in the red spectral domain. 

4.5. Periodicity peaks 

The histogram of the distribution of the variable sources 
versus their period is shown in Fig. The distribution 
of O-rich stars has a main peak at roughly Log(P) = 
2.80 and two lower peaks at about Log(P) = 1.95 and 
Log(P) = 2.30. C-rich stars have well defined peaks at 



about Log(P) — 2.45 and Log(P) = 3. Note that most 
of the sources with 2.53 < Log(P) < 2.59 have been ex- 
cluded because their period determination is suspected of 
aliasing. 

Fig. 4 in Cioni et al. lf2lM)l shows that O-rich AGB 
stars in the LMC have a strong broad peak between 
Log(P) = 1.5 and 2.0, and a lower peak at about 
Log(P) = 2.45, while C-rich AGB stars are homoge- 
neously distributed from Log(P) = 1.5 to 2.5. The first 
LMC peak of O-rich AGBs is probably a consequence 
of the selection criterium used to define the two sam- 
ples. The SMC sample is extracted from sources with an 
emission in the mid-IR, these stars are brighter and with 
longer periods than the stars in the LMC sample which 
were selected from the DCMC database. The difference 
in strength is probably due to the mean difference in the 
AGB type between the two Clouds. There are more car- 
bon stars in the SMC which is on average more metal 
poor than the LMC (Cioni & Habing EHE|) • This effect is 
also important for the sources populating relation C which 
are predominantly C-rich in the SMC and associated to 
the peak at about Log(P) — 2.45. That the majority of 
accepted Miras in the SMC are carbon stars was first re- 
cognized by Lloyd Evans et al. I|1988J1 as the result of a 
search for large-amplitude red variables in the Radcliffe 
Variable Star Field. The longer time range spanned by the 
MACHO observations allows to detect periodicities above 
Log(P) = 2.8 which was about the limit reached in the 
LMC study using EROS light-curves. The large number 
of sources detected in the SMC at these periods are likely 
to be multi-periodic. Their secondary period will be part 
of relation A, B or C, thus they will enhance the number 
of O rich stars between Log(P) = 1.5 and 2.5. In more 
detail we found that in our sample, among those with a 
red light-curve and q > 15, 24 sources have flag = 6. 
These sources have 2.0 < Log(P) < 2.3 which, taking into 
account the differences between the various period deter- 
minations and the mean, gives a period ratio from 1.5 to 
2. 

The same sample of LMC stars analysed by Cioni 
et al. I)2()()l|l was afterwards studied by Lebzelter et al. 
lISUU!?!) using the AGAPEROS datasets. Ninety per cent 
of their sources were classified into: regular, SR, irregular 
and other variability classes. Their analysis covers a time 
window of up to 900 days and is thus incomplete for the 
longest periods, but it is more complete for sources with an 
amplitude of variation below 0.1 mag. The authors found 
that the period distribution of regular variables has two 
peaks: one at Log(P) — 2.0 and the other at Log(P) = 2.5. 
The latter is probably formed by Mira stars. SR variables 
peak at about Log(P) = 1.8. The peak at the shortest 
periods is probably the same as the one found here. The 
peak of mostly Mira stars corresponds to the location of 
the peak of the C-rich SMC stars. In the Baade's win- 
dow, Glass & Schultheis ((212021 found that most of the 
late-type M giants have Log(P) < 2.0; those few with a 
secondary longer period have about Log(P) > 2.4. The 
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Fig. 5. Colour-magnitude diagrams (J—Ks, Ks) for vari- 
able stars classified as Mira (a), SRa (b) and SRb (c) ac- 
cording to their periodicity and amplitude of variation, 
(d) is the same diagram for variable sources which could 
not be classified in one of the three categories because of 
low S/N or an extremely irregular light -curve. The ver- 
tical dotted line indicates the separation between O-rich 
and C-rich stars as in Fig. 

histogram that can be derived from their Fig. 2 is very 
similar to that of O-rich stars in the SMC (Fig.^J;). 

4.6. (J — Ks, Ks) diagram 

In a similar fashion to the LMC sample (Fig. 7 in Cioni 
et al. I2(JU1|) . variable stars of a different type and varia- 
bility class are plotted in the colour-magnitude diagrams 
(CMDs) in Fig. |S| The distinction has been made as fol- 
lows: Miras have A x 100 > 110 independently of their 
type, and a regular light-curve with a well defined single 
periodicity. Similar sources but with a smaller amplitude 
of variation are SRa stars. SRbs are those stars which 
indicate the presence of a secondary period. Note that 
in this last category there is no distinction based on the 
amplitude of variation. As in the LMC, Mira stars are 
uniformly distributed in the regions where optically thin 
O-rich and C-rich, and also more optically thick stars, 
are located (Fig. SRa and SRb are almost equally 
distributed within the two branches of O-rich and C-rich 
optically thin AGB stars. Fig.[SJi shows the distribution of 
the sources that could not be classified in one of the three 
categories either because the light-curve has a too low 
S/N or it is extremely irregular. Most of these stars (180 
out of 194) have A > 10, thus they are variable stars and 
lie above the TRGB, K s = 12.62 (Cioni et al.ESSJ. It is 
possible that stars considerably brighter than Ks = 10 are 



supergiants and that fainter stars at about {J — Ks) = 0.4 
are foreground objects. Schultheis et al. (|2000(l found va- 
riable star candidates at (J—Ks) < 1. They conclude that 
it is probably a result of uncertainties in the determina- 
tion of extinction. However, they find a displacement of 
the Bulge stars detected as variables in both J and Ks to 
bluer colours than the sequence indicated by all detected 
ISO stars in the diagram (LW2-LW3, LW3). It would be 
interesting to know if some of these stars are supergiants 
and of which type. 

Finally we conclude that all the sources in the sample 
are variable stars because they have amplitudes above or 
equal to 0.1 mag either in the RED or BLUE wave band. A 
similar conclusion was reached by Alard et al. I|2001[) for a 
sample of sources in Baade's window detected at 7 and 15 
/im. The emission at the longest wavelength is a signature 
of mass-loss rate. Using their dependence between Mbol 
and the mass-loss rate and assuming that the mass-loss 
rate of highly evolved AGB stars does not depend very 
strongly on metallicity (van Loon 2000), we derive that 
the sources in the present sample lose mass at a rate from 
lO^MQyr- 1 (M BOL = -4) to l^M^yr" 1 (Mbol = 
— 7), but most of the objects have Mbol — —5 which 
corresponds to a mass-loss rate of about 10 _6 MQyr _1 . 
Note that a complete study of the dependence of mass- 
loss includes temperature, mass and outflow velocity as 
well. 



4.7. (K s ~ LW 10, LW10) diagram 

The 7/im LW2 filter does not include a large dust contri- 
bution for non-Mira sources and it is mostly a measure 
of the photospheric emission. On the other hand, for stars 
with a relatively thick circumstellar dust envelope, this 
band is strongly affected by the presence of dust (Alard et 
al. l2TMl Omont et al. H333)l . The 12^m LW10 filter, which 
is very similar to the IRAS 12/im band, includes the 9. 7/im 
silicate and the 11.3/im SiC dust features. Thus the colour 
(K — LW10) provides a strong indication of the amount 
of circumstellar dust for the entire sample. 

Fig. El shows the colour magnitude diagram (Ks — 
LW10, LW10) for the sources in the present sample, dis- 
tinguishing between Mira and SR variables. Except for a 
few SR all the sources with (Ks - LW10) > 2 are Mira 
variables (Fig. From the location of some K- and 
M-type supergiant (SG) stars in the ISO-MCMS one de- 
rives that they have (Ks - LW 10) < 1 and LW10 = 8. 
They follow a sequence that extends almost vertically to 
fainter magnitudes, but the nature of those sources with 
LW10 > 10 is not confirmed. Bright O-rich LPVs would 
have (K s - LW10) = 3 and LW10 = 7. These stars 
would have a late M subtype and not many of them are 
expected in the SMC. A quantitative statement will be 
possible once the study of the LMC LPVs detected by 
ISO is completed. A few SRb stars populate the fainter 
part (Fig. 03, c). There are many supergiant stars along 
the whole sequence that were not classified in any of the 
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Fig. 6. Colour-magnitude diagrams (K 5 — LW10, LW 10) Fig. 7. Colour-colour diagrams (Ks — LW2, LW2 
for the same groups of sources as in Fig. El LW10) for the same groups of sources as in Fig. 



three categories (Fig. It is thus not surprising that 
their light -curve is more difficult to classify because of a 
low S/N and irregularities. Thus, except for a few stars 
overlapping the region occupied by SR and Mira stars we 
arc able to distinguish sources with a poor or non-regular 
light curve. A more careful look at the near-IR colour- 
magnitude diagram (Fig.JBJ confirms that the unclassified 
stars with (J — Ks) — 1 or a bluer colour are experienc- 
ing the supergiant phase because AGB stars have overall 
(J—K s ) > 1. This is quite a good criterion to discriminate 
between supergiant and AGB stars, though the distinction 
is better achieved in the combined near- and mid- infrared 
diagram (Fig. EJ. 

Most carbon stars are located at about (Ks—LW10) = 
1 and 8 < LW10 < 11 (Fig. EJ-filled points); the sequence 
of dust obscured C-type stars extends to redder colours 
but the LW10 mag follows a shallower slope such as (Ks — 
LW10) = 4 at LW10 = 8. On average C-rich stars are 
redder than O-rich stars. 

Planetary nebulae and emission line objects, included 
in the ISO-MCMS catalogue, also have red colours (Ks — 
LW10 > 2) but they have fainter magnitudes with re- 
spect to Mira stars and are roughly located in a parallel 
sequence. Three of these objects are present in Fig. EJi. 
The object at LW 10 < 8.5 and K s - LW10 > 4 is a 
known carbon star while the object with the same LW10 
magnitude and a bluer colour is a known variable M giant, 
because of its red colour we erroneously classified it as a 
carbon star. 



4.8. (K s - LW2, LW2 - LWW) diagram 

Fig. E| shows the colour-colour diagram of the variables 
identified in this study The O-rich sources with (LW2 — 
LW10) > 0.4 are AGB stars or red supergiants. These 
AGB stars have a significant circumstellar extinction. This 
is also the colour of obscured O-rich and C-rich stars, and 
those with the highest amplitude of variability have (Ks — 
LW2) > 1.5. O-rich stars without a considerable mass loss 
have approximately (LW2 — LW10) — but can have < 
(Ks — LW2) < 1. There is not a remarkable distinction 
between SRa and SRb pulsators. At about (Ks — LW2) > 
3.5 one can find PN and emission line objects. 

It seems (Fig. 0i) that we may have missed the 
classification of a few Mira variables, O-rich and C- 
rich SR stars. At this stage it is worth going back to 
the specific light-curves to find out the reason of this 
missed classification and to perhaps identify some pe- 
culiar objects. Among the objects with Ks — LW2 > 
1.5 (Fig. J005113.6-731035, a known carbon star 
(RAW 658), has a flat light-curve, while the light-curve 
of J005136. 5-732016, identified as a PN, is irregular. 
Another PN is J005157. 9-731421 which may have a re- 
gular light-curve with an amplitude of about 0.8 mag. 
The periodogram of the other objects, J004849.4-731627 
(emission line star), J005212.9-730852 (variable M gi- 
ant) J005304.7-730409 and J010020.6-730648 (BMB-W 
29 carbon star) shows a very low amplitude with flare-like 
events. 
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Fig. 8. Histogram of the differences between the 2MASS 
and DCMC magnitudes of all the sources in the sample. 

4.9. Miscellaneous 

SRb, as have been defined in this work, do not constitute a 
separate class of variables. They are clearly multi-periodic 
objects and from their photospheric and dust colours they 
do occupy the same regions of SRa and Mira stars. Until 
a confirmation of the specific theoretical models (such as 
those developed by Winters et al. ( 1994J) to explain the 
multi- periodicity of the light-curve variation by coupling 
of stellar pulsation with the dust formation) or follow up 
radial velocity observations (to infer the binary nature of 
the sources) are available, the most appropriate class to 
which assign these stars is that of general SR variables. 
This is also strengthened by the fact that most of the ob- 
jects in the sample indicate the presence of more than one 
periodicity. A sophisticated observational and theoretical 
analysis of the multiple periods is required to further ad- 
dress this point. Because we did not discriminate on the 
amplitude of SRb stars some of them do have Mira-like 
amplitude but are likely to have more than one significant 
period. 

Most of the DENIS sources in the present sample were 
observed in December 1996 while 2MASS observations 
date August 1998. This difference in time of about two 
years may allow us to recognize variable stars by compa- 
ring their detection in the near-IR bands of DENIS and 
2MASS. Fig. [SJshows the distribution of these differences 
for all the sources in the sample. We checked that the 
separate histogram of sources with an assigned MACHO 
periodicity is the same as that of sources with Flag= 5 (a 
period was not assigned but the source is probably a vari- 
able). The histogram does not account for the systematic 
difference between the 2MASS and DCMC magnitudes as 



measured by Delmotte et al. I|2UU2|) . Selecting only the 
variable stars with the largest amplitude produces a si- 
milar histogram. The displacement of the peak from zero 
indicates that the sources are variable stars. Note that 
the difference is approximately the same in the J and Kg 
bands. 



5. Conclusions 

The stars selected from the ISO-MCMS with a confident 
DENIS/2MASS counterpart have been cross-identified 
with the MACHO database. The extracted light-curves 
have been analysed to derive the main period and am- 
plitude of pulsation. The bolometric magnitude has been 
calculated by integrating the spectral energy distribution 
from the / DENIS to the LTU10 ISO wave band. This, in 
combination with the correlation of Alard et al. ( 2001 ) be- 
tween Mbol and the mass-loss rate, indicates that most 
of the stars in the present sample lose 10 _6 MQ?/r~ 1 of 
matter. A distinction between Mira and SR is obtained 
by applying an amplitude and regularity selection crite- 
rion. The histogram of the amplitude of pulsation of O- 
rich and C-rich stars has a similar distribution, but on 
average C-rich stars have a larger amplitude. This may 
indicate that either most of the C-rich stars in the SMC 
are of Mira type, or that the metallicity affects the ampli- 
tude in such a way that in a lower metallicity environment 
the amplitude of pulsation is larger. This effect cannot yet 
be checked for in other metal-poor galaxies in the Local 
Group, because despite the fact that many AGB stars have 
been discovered there is not enough information on their 
variability and type. 

The most obscured stars have large (Ks — LW2) and 
(LW2 — LW10) colours. This indicates sources with a 
thicker circumstellar envelope and allows to easily distin- 
guish AGB stars from supergiant stars and other classes of 
objects such as PN and emission line sources. The differen- 
ce between the 2MASS and DCMC magnitudes suggests 
that most of the sources in the sample are variables. 
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Appendix A: The period-luminosity relation of 
Mira stars 

In order to improve our understanding of the identification 
of Mira-type variables we decided to re-analyse the LPVs 
studied by Feast et al. 119891 These are the Mira stars 
that define the period-luminosity relation in the LMC. 
We could find the B1950 coordinates of 42 out of the 55 
sources with a "good" or "fair" light-curve as classified 
by Feast et al. I|1989|l . We could not find in the literature 
the coordinates of the objects from the sample of Glass & 
Lloyd Evans (1981) and for three objects from the list of 
Glass & Reid fiSEty or Reid et al. lTI$ggjl . However, there 
still remain enough sources for our purpose. 

We searched in the MACHO database for the light- 
curve of the closest LPV to these 42 sources. The iden- 
tification was easy in most of the cases because of the 
highly regular and large amplitude light-curve typical of 
a Mira star. Five sources (W158, GR28, GR4, RGC91, 
RGC89) 3 are unfortunately outside the fields observed by 
MACHO. For one source the MACHO web site gives a 
persistent initialization error (W94) and for four addi- 
tional sources (GR13, GR21, RGC37, RGC53) we could 
not find a nearby counterpart (a first possible candidate 
is further than 8"). Only one source (W30) has two pos- 
sible MACHO counterparts and due to the low accuracy 
of the initial coordinates compared to the MACHO coor- 
dinates we decided to remove this source from the discus- 
sion. Finally we could successfully extract the light-curves 
of 31 sources. 

We determined the main periodicity as described in 
Sect. 3, of the sources with a MACHO-red light-curve. We 
obtained periods that are at most 20 days different from 
those derived by Feast et al. I|1989)) except for 4 sources 
which result to have a period coincident with one of the 
extremes of the searching window (GR2, GR3, GR18 and 
GR29, they have been removed from the following analy- 
sis) and source GR17 for which we derive a period of 51 
days shorter. The period of source RGC20 falls in the re- 
gion of aliases (Sec. 3), thus the source has been removed 
from the following discussion. The quality parameter of 
the 26 sources with a convincing period determination is 
always above 100 and the difference between the assigned 
period and the mean of all the determinations with q > 15 
is of the order of 1 day for most of the sources but always 
below 10 days. These numbers may be different in the case 
of MACHO-blue light-curves. 

Tablc lA~Tl lists the name of the sources (Column 1), the 
coordinates at the epoch J2000 precessed from the original 
coordinates (Columns 2 & 3), the MACHO identification 
code (Columns 4-, 5 & 6), the distance in arcsec between 
the Feast et al. (1989) source and the MACHO counter- 
part (Column 7), the period, amplitude and the value of 
the quality parameter for the red light-curve ( Columns 8, 
9 & 10). 

3 The designation of the sources follows that of Feast et al. 
IT3551 W=Wood et al. 115551 GR=Glass & Reid YFMA and we 
introduce RGC=Reid et al. 115551 
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Fig. B.l. The bolometric correction BC/ as a function of 
the I — J colour for the O-rich (open circles) and C-rich 
(filled circles) variable stars in the sample. The continuous 
line is the relation derived by Alvarez et al. (2000) for cool 
O-rich stars. 

We conclude that by selecting, from our sample, large 
amplitude sources with a single periodicity and q > 100 
we are identifying Mira variables of the same type as those 
studied by Feast et al. 1989 They define a tight PLR and 
the small scatter in Fig. is probably due to the availa- 
bility of only one epoch bolometric magnitude. Note that 
with this selection criterion we also find sources at very 
bright magnitudes (Mbol = ~~ 7) which deviate from the 
extrapolation of the PLR defined only for sources with 
P < 420 days. We estimate from Fig. that these stars 
have M > 5Mq and as previously noted by Feast et al. 
1989 the PLR breaks at high masses. These stars are likely 
to be "over-luminous" as a result of Hot Bottom Burning 
(Zijlstra et al. GSHSl Whitelock et al.[2003). However, be- 
cause of their very high luminosity they could also be su- 
pergiants. 

Appendix B: Bolometric correction(s) 

B.l. BQ versus I - J 

Fig. IB . ll compares the bolometric correction obtained from 
the integration under the SED (Sect. 3.1) and the relation 
derived by Alvarez et al. |2000) for cool O-rich stars. The 
agreement is very good for the O-rich stars of our sample. 
This confirms that the relation is valid for the SMC as well 
as for the LMC and the Bulge. However this relation over- 
estimates the correction for C-rich variables in the SMC. 
On the other hand we should check that the integration 
under the SED does not underestimates the bolometric 
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correction (if the stars have a circumstellar envelope and 
dust extinction the correction for energy re-radiated by 
dust becomes important). Our measurements have been 
corrected for interstellar extinction while the relation has 
not. The agreement for O-rich stars might be a conse- 
quence of the systematic difference between the DENIS 
filters, used in this work, and the broad pass-bands by 
Bessels & Brett (1988| and Bessel (|1990fl used to derive 
the relation. The disagreement for C-rich stars might be 
a consequence of the SiC dust feature which causes the 
spline not to follow the continuum of the SED. However 
the data of Alvarez et al. l|2(JU(Jfl include only 6 C-rich 
stars and a total of four stars belonging to the LMC or 
the SMC. In addition the spectral range covered by their 
measurements is limited to 2.5/im. The zero-point of the 
bolomctric scale is the same in both studies. 

B.2. BC K versus J - K 

Fig. IB. 21 compares the bolometric correction obtained as 
in Sect. 3.1 with the relation derived by Montegriffo et al. 
ifrM^ for Population II giants of low ([Fe/H]< -1.0) and 
"high" metallicity ([Fe/H]> -1.0). They used M BOLj3 = 
4.75. Their relation was defined up to (J — K) = 0.8 but 
being the metal-poor relation very close to the metal-rich 
relation at redder colours we consider the latter as a pos- 



sible extrapolation. This results in a fairly good fit for 
the O-rich variables while C-rich variables lie systemat- 
ically below the relation. This is a similar effect as that 
observed for BCj versus I — J (Fie. IB.ljl . The giants in 
the sample of Montegriffo et al. (1998) are members of 
galactic globular clusters. C-rich stars are rare members 
of these systems. The spectral range covered by the ob- 
servations is limited to 3.5^m (L band) and they derive 
the bolometric magnitude by integrating the SED which is 
previously re-constructed by means of a set of Plank func- 
tions at a given wavelength. These effects may influence 
the correction to be applied to the flux measured only in 
the K band. Besides the 2.2/im filters used in both studies 
are not strictly the same. We conclude that we cannot at 
present compare the correction for C-rich stars while O- 
rich stars are qualitatively in agreement. It would be useful 
to compare stars of other galaxies in the Local Group at 
a different metallicity. 
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